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In the framework of the Chandrasekhar mass white dwarf model for Type la supernovae, various 
stages of the explosion are described in terms of the burning regimes of the thermonuclear flame 
front. In the early flamelet regime following the "smoldering" phase prior to the explosion, the 
flame is sufficiently thin and fast to remain laminar on small scales. As the white dwarf density 
declines, the thermal flame structure becomes subject to penetration by turbulent eddies, and it 
enters the "distributed burning" regime. A specific control parameter for this transition is proposed. 
Furthermore, we outline an argument for the coincidence of the transition between burning regimes 
with the onset of a deflagration-detonation-transition (DDT) in the late phase of the explosion. 



I. INTRODUCTION 



Theoretical efforts to understand the physics of Type la supernovae (SN la's) based on the model of exploding 
Chandrasekhar mass white dwarfs have revealed a fascinating degree of complexity. While the basic model appears 
simple - a carbon and oxygen white dwarf that, by some means, reaches critical conditions for thermonuclear runaway 
inevitably burns a sizeable fraction of its mass to 56 Ni, releasing roughly the amount of energy needed to power the 
explosion (~ 10 51 ergs) - the hydrodynamics of the explosion is far from being fully understood. As is well known, 
, the observed abundance of intermediate mass elements in the ejecta rules out the possibility that the star is disrupted 
by a single, prompt detonation (a shock-driven combustion wave), which would burn almost the entire star to nuclear 
s ! \ statistical equilibrium [Q . Consequently, any successful model must involve an initial - conductively propagating and 
thus subsonic - deflagration, or "flame", phase that allows the star to pre-expand. Turbulence driven by buoyancy 
of the burnt material with respect to the unburnt background increases the total burning rate, first by wrinkling the 
flame surface and later by directly mixing hot ashes with cold fuel. Only allowed to occur below a critical density 
of ~ 10 7 g cm -3 , the process of direct turbulent mixing might be responsible for triggering a detonation in the late 
, stage of the explosion or during the re-collapse of the partly burned white dwarf H||. 

The complicated interaction of thermonuclear burning, thermal conduction, turbulent mixing, and buoyancy forces 
on scales covering up to 12 orders of magnitude makes a full numerical solution of the governing equations virtually 



I ■ impossible. Multidimensional simulations on the scale of the white dwarf radius suffer from insufficient resolution of 
small scale turbulence and of the thermonuclear flame front, and thus (implicitly or explicitly) involve subgrid modeling 
+3 . to some extent, rendering the interpretation of their results intrinsically difficult. One of the central assumptions made 
in previous calculations is the notion of a thin flame surface separating fuel (in our case, carbon and oxygen) and 
ashes (nickel), and propagating into the fuel at a speed that only depends on its density and composition. Under 
certain conditions, the laminar flame structure remains unperturbed by turbulent eddies on small scales, while at 
the same time turbulence wrinkles or even fragments the flame surface on larger scales. This so-called "corrugated 
flamelet regime" Q occurs early during the explosion, when the flame is fast and thin and turbulence is weak. Key 
parameters for large-scale simulations in the flamelet regime are the boundaries of validity of the flamelet assumption, 
the small scale cut-off for flame surface perturbations, and the turbulent flame speed (the global propagation rate of 
the turbulent flame brush) as a function of the laminar flame speed and turbulence intensity. Section (fl| ) gives an 
overview of thermonuclear burning in the flamelet regime. 

It has recently become obvious that one cannot ignore the effects of turbulence on the laminar flame structure once 
the star has expanded to densities below ~ 10 7 g cm" 3 ||. As soon as turbulent eddies on the scale of the flame 
thickness become sufficiently strong to disrupt the thermal flame structure and mix fresh fuel into the hot burning 
products, the flame enters the "distributed burning regime" ||. It has been noted ||] that this process generates 
favorable conditions for a deflagration-detonation-transition (DDT) that is thought to occur in many SN la events 
. Clearly, an explicit control parameter that describes the transition from flamelet to distributed burning would be 
useful; one possible choice, the ratio of turbulent and conductive diffusivities, will be discussed in section (|Tv|). 

The following three sections summarize the general features of the various stages of a Chandrasekhar mass SN la 
explosion, described in terms of their respective burning regimes. It should be kept in mind that the relations sketched 
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below are chiefly based on dimensional and scaling arguments. Wherever feasible, accurate numerical simulations will 
have to be carried out in order to verify or discard these arguments and to specify numerical coefficients. 



II. PHASE 1: SMOLDERING AND IGNITION 



Very little is known about the time between the onset of the thermonuclear runaway and the formation of the flame 
itself. At a central temperature of T 2 x 10 8 K, neutrino cooling fails to keep up with nuclear heating, and the core 
region begins to "smolder" [raj. During the following ~ 1000 years, the core experiences internally heated convection 
with progressively smaller turnover time scales r c ~ v7 , where Vbv is the Brunt- Vaisala frequency. Simultaneously, 
the typical time scale for thermonuclear burning, 7b, drops even faster as a result of the rising core temperature 
and the steep temperature dependence of the nuclear reaction rates. At T w 6 x 10 8 K, both time scales become 
comparable, indicating that convective plumes burn at the same rate as they circulate Experimental or numerical 
data describing this regime of strong reactive convection is not available. 

At T rs 1.5 x 10 9 K, Tb becomes extremely small compared with t c , and carbon and oxygen virtually burn in place. 
A new equilibrium between energy generation and transport is found on much smaller length scales, I rs 10~ 4 cm, 
where thermal conduction by degenerate electrons balances nuclear energy input The flame is born. 

The evolution of the runaway immediately prior to ignition of the flame is crucial for determining its initial location 
and shape. Using a simple toy model, Garcia & Woosley [p"T| found that under certain conditions, burning bubbles 
subject to buoyancy and drag forces can rise appreciable distances before flame formation, suggesting the possibility 
of off-center ignition. As a consequence, more material burns at lower densities, thus producing higher amounts 
of intermediate mass elements than a centrally ignited explosion. A parameter study demonstrated the significant 
influence of the location and number of initially ignited spots on the final explosion energetics and nucleosynthesis 

& 

The topology of the initial flame surface is directly linked to the spatial and spectral state of temperature fluctua- 
tions due to the strongly coupled dynamics of three-dimensional convection in a stratified medium, microscopic heat 
transport, viscous dissipation, and nuclear burning. Presently, there is no obvious way to neglect or parameterize any 
of these processes in numerical simulations without risking qualitatively wrong results. 



III. PHASE 2: THE FLAMELET REGIME 



The general features of laminar thermonuclear carbon and oxygen flames at high to intermediate densities were 
described in detail by Timmes & Woosley |10[j . For our purposes, we need to know the laminar flame speed uq w 
10 7 . . . 10 4 cm s _1 for p rs 10 9 . . . 10 7 g cm -3 , the flame thickness ith — 10~ 4 ... 1 cm (defined here as the width of the 
thermal pre-heating layer ahead of the much thinner reaction front), and the density contrast between burned and 
unburned material p = Ap/p = 0.2 . . .0.5 (all values quoted here assume a composition of Ac = Ao = 0.5). The 
thermal expansion parameter p reflects the partial lifting of electron degeneracy in the burning products, responsible 
for the transformation of pre-ignition convection into a genuine Rayleigh- Taylor (RT) problem after formation of the 
flame. 

In the standard picture, buoyant bubbles of ashes rising through the fuel create turbulent velocity fluctuations v! 
on the scale of their diameter d, which cascade down to smaller scales. By thermal expansion and laminar flame 
propagation, bubbles grow to d ss 10 6 . . . 10 7 cm during the first 10 _1 s of the explosion. Subject to the gravitational 
acceleration g ss 10 9 cm s -2 , their terminal rise velocity is 

v r ~ 0A(pgd) 1/2 w 10 7 cm s" 1 . (1) 

Hence, Uq/u' ~ uo/v r <C 1 and lth/d 1. These are two necessary conditions for the flamelet regime; a third 



condition relating the turbulent diffusivity on small scales to the thermal conductivity will be discussed in (IV). Since 
the propagation of the turbulent flame brush is dominated by the velocity of the largest turbulent eddies, the turbulent 
flame speed «t becomes independent of uq, and 

UT „ u > if ^ _> o (2) 
w 

follows from simple scaling arguments. This behavior is strongly supported by experiments Jl3| . In our context, 
a short time after flame ignition the total burning rate of the turbulent flame brush essentially decouples from the 
microscopic physics of nuclear burning and heat transport, and from there on depends only on the hydrodynamics of 
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buoyancy-driven turbulence. Hydrodynamical simulations of the flamelet phase in SN la explosions make use of this 
argument (see below). 

Given the smoothness of the flame on scales Zth and the existence of large perturbations on scales d Zth, there 
must be an intermediate scale l cu toS corresponding to the transition between both regimes. In turbulent chemical 
combustion, where turbulent velocities generally scale according to Kolmogorov scaling, u'{l) ~ I 1 / 3 , this scale is 
known as the "Gibson length" One can find an estimate for Z C utoff by assuming that the flame surface is only 
affected by eddies that turn over at least once during their flame crossing time: r ddy = fcutoff/w' (Cutoff) = fcutoff/iio, 
and therefore u'(Z C utoff) = uq. For v! = u'(d) and Kolmogorov scaling, the cutoff scale for flame surface perturbations 
is /cutoff = d(u Q /u') 3 . 

The same argument holds true for flames in SN la explosions, except that here the turbulence that deforms the 
surface is created by buoyancy of the hot burning products and is not, as in most terrestrial experiments, generated 
by a grid. As a result of the parallel cascades of kinetic and potential energy (due to the presence of both velocity 
and density perturbations), the resulting velocity spectrum differs from Kolmogorov scaling. Rather, buoyancy-driven 
turbulence conforms to Bolgiano-Obukhov scaling JTif , 

u'{l) ~ l 3/b . (3) 
The cutoff scale relevant for the analysis of the flamelet regime in SN la explosions is thus 

W-rf(^) 573 . (4) 

Below the cutoff scale, the flame surface is smooth. Buoyancy does not couple to the turbulent kinetic energy cascade 
on these scales due to the absence of density fluctuations. Hence, the velocity spectrum turns over to Kolmogorov 
scaling at i C utoff > continuing down the viscous microscale. 

We can use this piece of information to specify the requirements for hydrodynamical simulations. As long as Z C utoff 
is unresolved, the propagation velocity of the numerical "flame brush", m(A), is not equal to w ; m fact, @ suggests 
that it becomes independent of Uq in the case of strong turbulence. One possible alternative is to employ the rough 
approximation ut(A) ~ u'(A), where the turbulent velocity on the grid scale A can be extracted from a subgrid 
model for the turbulent cascade Jl^]. However, with typical values for u' ~ 10 7 cm s _1 and d ~ 10 7 cm, Cutoff is not 
microscopically small at high densities (~ 10 4 cm for p w 10 8 g cm -3 ), and full resolution of this scale is within reach 
of future simulations. In this case, the use of uq for the flame propagation rate is justified. The unphysical situation 
associated with not resolving the viscous cutoff, of course, remains. 

As a final remark on the validity of Gibson scaling, ? cu t ff ~ (uq/u 1 ) 3 in Kolmogorov turbulence, the turnover 
to smoothness at this scale has not yet been confirmed experimentally or by direct simulations. This fact has 
been attributed to flame-turbulence interactions by thermal expansion that are absent in the passive-surface picture 
employed above. In white dwarf explosions, partial degeneracy of the burning products results in a much lower degree 
of thermal expansion than in chemical flames, where fi is generally of order unity. Therefore, the passive-surface 
approximation is, compared to terrestrial flames, much less restrictive. For purely passive flames, Gibson scaling was 
clearly demonstrated in a simple discrete model for turbulent flamelet combustion Jl6[ | that achieved spatial resolution 
of the flame surface over four orders of magnitude. 

Assuming that buoyancy effects dominate the turbulent flow that advects the flame, the passive-surface framework 
obviously neglects the additional stirring caused by thermal expansion within the flame brush itself, accelerating 
the burnt material in random directions. Both the spectrum and cutoff scale may be affected by "active" turbulent 
combustion ||. Although the small expansion coefficient /i indicates that the effect is weak compared to chemical 
flames, a quantitative answer is still missing. 



IV. PHASE 3: THE DISTRIBUTED BURNING REGIME AND POSSIBLE 
DEFLAGRATION-DETONATION-TRANSITION (DDT) 

As the density of the white dwarf material declines and the laminar flamelets become slower and thicker, a point is 
reached where turbulence significantly alters the thermal flame structure. This marks the end of the flamelet regime 
and the beginning of the distributed burning, or distributed reaction zone, regime. In order to find the critical density 
for the transition between both regimes, we need to formulate a specific criterion for flamelet breakdown. 

In the combustion literature, the Karlovitz number Ka = ith/^viso where l v i sc is the viscous cutoff scale, is commonly 
used to characterize the transition jlTj . The flamelet regime corresponds to Ka < 1, while Ka > 1 implies the existence 
of turbulent eddies smaller than the flame thickness, which is often interpreted as the onset of distributed burning. An 
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alternative definition of the Karlovitz number compares the dissipation and reaction time scales, Ka = Tb/T v ; sc [fL8[ . 
As is readily seen by inserting Tb = r t h = lth/ K anc ^ T visc = ^visc/ 1 ^ where k is the thermal diffusion coefficient and v is 
the kinematic viscosity, these definitions are only equivalent for Prandtl numbers Pr — v] k near unity. In contrast to 
chemical flames where Pr « 1 is usually a reasonable approximation, the Prandtl number in the degenerate medium 
of a Chandrasekhar mass white dwarf is much smaller, Pr ss I0~ 5 . . . lO -1 |19|. The length and time scale criteria 



diverge in this case, implying that the conditions hh/lvisc > 1 and Tb/r V i SC < I can coexist. If eddies smaller than 
the flame thickness exist, but they are completely burned before they can turn over, the Karlovitz criterion lacks an 
obvious interpretation. 

A transition criterion that is independent of Pr was proposed by Niemeyer & Kerstein || . On phenomenological 
grounds, one can argue that if the ratio of turbulent diffusivity, K e (l) ~ lu'(l), and microscopic heat diffusivity 
k exceeds unity on the scale l t h or below, the flamelet regime breaks down. Since the turbulent diffusivity is an 
increasing function of scale, it is sufficient to evaluate k c on the scale lth- The diffusivity criterion can then seen to 
be equivalent to 

" > 1 (5) 



Teddy (^th 

for flamelet breakdown. Dividing both time scales by Z t h, one finds that it coincides with u'(l t h) > uq and hence with 

^cutoff < 'th ■ (6) 

Inserting the results of Timmes & Woosley for uq and lth as functions of density, and using a typical turbulence 
velocity u'(I0 6 cm) ~ 10 7 cm s _1 , the transition from flamelet to distributed burning was shown to occur at a density 
of p « I0 7 g cm -3 H. Intriguingly, one-dimensional SN la models that artificially invoke the onset of a deflagration- 
detonation-transition (DDT) after a slow initial flame phase achieve best agreement with observations if the DDT 
occurs very close to this density [@]2(]]. Lacking a physical description of the DDT itself, this fine-tuning of the 
transition density is unnatural. 

If, however, the conditions for DDT become more favorable in the distributed burning regime than in the flamelet 
regime, the initiation of a detonation could be naturally related to flamelet breakdown. It was proposed jj| that the 
local quenching of flamelets induced by turbulent mixing of fuel and ashes can lead to the formation of macroscopic 
"smoldering" regions, i.e. material that still burns, but on a much longer time scale than a laminar flame under 
the same conditions. If only one such region re-ignites in a nearly isothermal state, brought about by continuous 
turbulent stirring, it can evoke a detonation. In order to re-ignite before the material has expanded substantially, 
Tb in the smoldering region must be smaller than the dynamical time scale for expansion, Td ~ 0.1 s. On the other- 
hand, turbulent mixing must be sufficiently fast to homogenize the fluid over a length scale given by the critical 
mass for detonation, Z cr ;t ~ (m C rit/ p) 1 ^ 3 '- if the the near- isothermal region is too small at the time of ignition, the 
ensuing pressure wave fails to trigger a self-sustaining detonation. The critical mass (or length) depends sensitively on 
composition, density, and boundary conditions [^|,^) and is not well known. However, under reasonable assumptions 
the ordering of time scales T e ddy('crit) < Tb < Td is fulfilled ||. 

Even if the probability for a single DDT were very small, this may be compensated by the (potentially) large 
number of critical masses in the flame brush. Detailed knowledge of m cl -it under various conditions, determined by 
direct numerical simulations of turbulent smoldering fluid regions followed through re-ignition, may eventually enable 
us to estimate the global probability of a DDT in the distributed burning regime. 



V. SUMMARY AND CONCLUSIONS 



The main goal of this essay is to classify the phases of the thermonuclear explosion of a Chandrasekhar mass white 
dwarf in terms of the burning regimes experienced by the turbulent combustion front. By comparing the various 
length and time scales characterizing the flame, three major regimes can be identified: the "smoldering" regime prior 
to the dynamic part of the explosion, followed by the flamelet regime where the flame structure is microscopically 
laminar but corrugated on large scales, and finally the distributed burning regime, occurring when turbulent transport 
on the scale of the flame thickness begins to dominate over microscopic heat diffusion. 

Of these three, the flamelet regime is probably the most studied and best understood. As long as the "passive 
surface" and "thin flame" assumptions hold (as specified in the main text), the turbulent flame speed decouples from 
the physics of nuclear burning and microscopic heat transport. Instead, it scales with the velocity of turbulent eddies, 
which can be estimated using the rise velocity of hot bubbles and the Bolgiano-Obukhov spectrum of buoyancy-driven 
turbulence. The full problem involving feedback of thermal expansion on the turbulence ("active combustion"), time- 
dependent gravitational acceleration, freezing-out of the large scale eddies due to expansion of the star, and parallel 
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cascading of kinetic and potential energy continues to be accessible only to extensive numerical simulations; however, 
the orders of magnitude of the governing parameters are probably already known to reasonable accuracy. 

The smoldering and ignition phase, on the other hand, is only poorly understood. Regarding the evolution of the 
explosion itself, the most important information we would like to extract is the initial shape and location of the 
flame front. Likewise, the question whether the flame is turbulent right from the beginning or makes a transition to 
turbulence from an initially smooth state is still unanswered. Detailed numerical experiments that fully resolve the 
microphysical transport and burning processes are required to address this subject. 

Most important from the point of view of SN la modeling is the transition from flamelet to distributed burning as 
the density drops below 10 7 g cm -3 , as this may present a physical mechanism for the onset of a delayed detonation. 
Confirmation of this mechanism would solve one of the outstanding fine-tuning problems of one-dimensional SN la 
models that, if the turbulent flame speed and the density for the deflagration-detonation-transition (DDT) are suitably 
chosen, agree well with observed SN la spectra and lightcurves. Further understanding of the physics of re-ignition is 
also needed to decide whether the DDT occurs during the first expansion phase or after one or several pulsations. Since 
the proposed mechanism for DDT involves local flame quenching and re-ignition in a homogenized state, the set-up 
for hydrodynamical simulations is very similar to the initial ignition problem, only at a higher degree of background 
turbulence and lower density. Clearly, the dynamics of a turbulent, slightly sub-critical medium warrants closer study. 

Given the uncertainties in the ignition and DDT processes, one can easily find possible explanations that account 
for inhomogeneities among the observed SN la events. However, lacking a detailed understanding of the underlying 
physics, it is difficult to evaluate how sensibly the global outcome of the explosion depends on minor variations of the 
initial conditions or the transition criteria for burning regimes. 
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